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Abstract 

Biological systems are among the most challenging subjects for theoretical physicists, 
as well as experimentalists or simulationists. Physical principles should have been both 
constraints and guide-lines for the evolution of living systems over billions of years. One 
generally aims at clarifying what physical principles, possibly new ones, are behind the 
phenomena of biological interest and at understanding how they work within the entire 
biological world. In the present talk we describe an example of such an effort. 

Since the discovery of 'superprecipitation' by Szent-Gyorgyi 's group in 1940's, it has 
been a long puzzle how an assemblage of actin filaments with random orientation can 
contract in the presence of the two-headed myosin molecules undergoing actin-activated 
ATP-hydrolysis reaction. It is widely accepted that during the contraction the two-headed 
myosin mediates the relative sliding of two actin filaments whose polarity directions are not 
parallel but rather anti-parallel. But this fact solely does not account for the shortening. 
We propose a dynamical model which, upon numerical simulation, exhibits the shortening 
of an bundle of the actin filaments which are initially dirstributed randomly both in space 
along a line and in polarity direction. In the course of shortening several clusters of 
actins appears along the bundle. The model also shows the sorting of the actin filaments 
according to their polarity in the late stage. These findings are in accordance with the 
recent experiment by Takiguchi. 
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1 Introduction 



There is much interest in biological systems from a physicists point of view by several 
reasons. First, by looking at the diversity and hierarchy of biological phenomena and at 
the billions of years of their evolution, it is a challenge to unveil universal phenomena 
or universal origins of those systems based on physical principles. For instance, ATP 
(adenosine triphosphate) is often described in biology textbook as the energy source [energy 
donner] , as the substrate of transferase [phosphate donner] , and as the substrate of allosteric 
enzyme [regulation factor] . We are, however, tempted to search more unified view of the role 
of ATP since it should have appeared on the earth initially bearing a single role. Secondly, 
the biological systems are the very representatives of complex systems. By regarding them 
as systems of active elements we are inspired with many physical ideas and models. 

If we view the subjects of biology which have become also the subjects of physicists, 
we find that there are some frameworks elaborately introduced so that physicists can de- 
velop their idea upon it. Protein folding is studied based upon the Anfinsen's dogma that 
(most) natural proteins rest in their equilibrium folded states. Neural network, despite the 
prohibiting complexity in reality, is studied based upon mathematical realization of Hebb's 
hypothesis. Fluctuating membrane is studied on the basis of elasticity theory including 
entropic or Helfrich interaction, and molecular evolution is studied as stochastic process. 
Morphogenesis and pattern formation have been studied in the framework of bifurcation 
theory, etc. Protein dynamics can be one of the near future target of physicists, being 
stimulated by the recent development and need of nanoscale handling of soft materials, 
though the theoretical framework for it is not yet established. 

It is said that the biological processes which appear to be purely physical phenomenon, 
such as symmetry breakdown or instability are even sometimes coded explicitly on the DNA. 
It should be, however, still meaningful to ask "how did such a biological process happen 
to appear and become incorporated into evolutionary process?" Upon the appearance of 
a new biological function, it should have been quite primitive and unsophisticated, which 
works in barely efficient way or in a poorly organized way. Exploring the mechanism of 
such primitive functions should then be a subject of physics of biological interest as well 
|]1[]. As such an example, we present in this paper how the system of random assemblage of 
myosin and actin filaments (both being the constituent proteins of our muscle) can exhibit 
stochastic contraction phenomena, which is recently studied experimentally in detail [Q]. 
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Figure 1: 
(a) A head of the two- 
headed myosin ('seed- 
leaf') translocates an 
actin filament (thin ar- 
rowed line) in the for- 
ward direction indi- 
cated by an open thick 
arrow. (b) and (c) 
The action of a myosin 
onto two filaments. 




Figure 2: To'p: Uniaxial assem- 
blage, or bundle, of actin fila- 
ments, and the myosin molecules 
bound to them. Bottom: As the 
myosins translocate the filaments, 
overall shortening of the bundle 
occurs. 



Figure 3: The action of 
myosin can contributes to 
either (a) shortening or (b) 
to elongation of the actin 
bundle. In the case (b) it 
can happen that the elon- 
gation is interrupted, as ar- 
gued in the text. 



2 System of many myosins and actin filaments: 
A paradox 

We are considering the system consisting of myosins and actin filaments. Each myosin 
has two globular heads (shown by the symbols like seed-leaves in Fig. |l]), and each actin 
filaments has its own polarity direction (indicated by the thin arrow lines in Fig. |lD. If a 
globular head of a myosin is within the reach of an actin filament in the presence of ATP, 
the globular head consumes the hydrolysis energy of ATP to drive the actin filament to the 
forward direction (indicated by the open thick arrow in Fig. |l](a)). The net relative motion 
of actin filaments is brought by the action of a myosin only when the myosin is bound to 
non-parallel pair of filaments (Fig. |l](c)), but not to parallel pair (Fig. |l](b)). 

We focus here on the recent experiment by Takiguchi 0, He prepared a bundle of 
many actin filaments which are assembled uniaxially but randomly with respect both to 
their position and to their polarity direction. It has been demonstrated [H that the bundle 
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of actin filaments undergoes longitudinal shortening in the presence of myosin and ATP 
molecules. We will describe the experimental procedure of p[ in more detail. First a 
long and thick bundle of many actin filaments is prepared in methyl cellulose aqueous 
solution. To this bundle the two-headed myosins (so called heavy meromyosin, or HMM) 
and an abundance of ATP molecules are added. The bundle then starts to contract slowly 
in length, while it thickens so as to conserve its volume (Fig. |^). This shortening often 
occurs in a wiggling way. After the bundle has shortened appreciably, several needle-like 
subbundles appear from the main bundle. It has been shown that in these subbundles the 
polarity of the actin filaments is not random but is oriented outward with respect to the 
original bundle. 

Experiment like this dates back to late 40's, when Szent-Gyorgyi's group discovered so 
called superprecipitation, the phenomenon that a three-dimensional random assembly of 
actin filaments and myosins shrinks dramatically after the addition of a certain amount of 
ATP 0. Such experiment has been recently also done and refined Takiguchi's setup 
may be regarded as a more idealized one to see how the contraction occurs. Such an ideal 
random distribution is realized only in in vitro experiment, but the situations more or less 
like this have been found in nature such as in the contractile rings that appear during the 
mitotic period of cell division cycle |0 or in stress fibers observed in the locomoting cells 
during their contraction period 0. The experiments mentioned above, therefore, could 
be regarded as a simulation of in vivo systems or, at least, as a hypothetical simulation 
of evolutionally primitive stages of muscle contraction or cell motility. The question how 
this primitive system undergoes shortening has, however, not been studied for a long time 
since the discovery by Szent-Gyorgyi's group. It is because the highly organized structure 
of muscle was found [§] shortly after the former discovery, and the main stream of muscle 
study has been focused towards a dynamics of single globular head of myosin and its 
regulation mechanism 0. 

It is Hatano who seriously questioned how the actin bundle can shorten in the prim- 



itive situation like in Fig. ^, and he came across the following paradox |10]: When the 
sliding of the oppositely oriented actin filaments starts from the state shown in Fig. ^(a), 
the overlap between the two filaments would increase, leading to the shortening of the 
bundle. On the other hand, when the sliding of the filaments starts from the state shown 
in Fig. Wb), the action of the myosin would decrease the overlap between the filaments. 
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leading to the elongation of the bundle. Since the both situations should occur equally 
likely in a bundle, there would be no net shortening at all. In fact so-called bipolar kinesin, 
the other motor protein closely related to myosin, is discovered to appear during the cell 
division process and this protein is thought to act to separate the two spindle-poles by the 
mechanism shown in Fig. ^(b) |^. We would note that the above paradox cannot be lifted 



by considering the effect of simultaneous action of many myosins to an actin filament, as 
it occurs experimentally as far as we assume the continuous action of myosin molecules on 
the actin filaments, while such model could predict the undulational instability of filament 
density 



3 Simple model and simulation 

Our simple idea to resolve Hatano's paradox is to take into account the finite distance, 
say b, by which a globular head of myosin can continuously drive a single actin filament 
(Fig. 1^). The limitation of this distance may come from the dynamic fluctuation of the 
myosin heads as well as by the fluctuation of the lateral arrangement of actin filaments 
within a bundle. Our reasoning for the shortening is as follows: If the myosin acts in the 
situation of Fig. |^(a) each globular head can translocate the respective filament fully by 
the distance b on the average (hereafter we assume that b is sufficiently smaller than the 
length of each filament, which we denote by i), while in the situation of Fig. ^(b) the 
translocation of actions by myosin can be interrupted when one of its two heads meets 
with the rear end of an actin filament. The elongation of the bundle to which the filaments 
shown in Fig. ^(b) belong is, therefore, less extensive than the shrinking of the bundle 
to which the filaments shown in Fig. |^(a) belong. The interruption of the elongation will 
occur by the probability proportional to b/i in the approximation up to the lowest order 
of 6/£ <^ 1. The net shrinkage per each action of myosin will then be roughly scaled by 
~ 6 • 6/^ if a single myosin acts to the pair of filaments. Actually the shortening by this 
mechanism should be still less efficient due to the presence of other myosins interacting 
with those actin filaments. We believe, however, that the basic mechanism of shortening 
may be captured by the present simple model. 

We performed a numerical simulation based on the simple idea described in Fig. ^. The 
algorithm of the simulation is as follows: First we distribute randomly N actin filaments of 
the length £ over an interval —^<x<^ along the a;— axis. We choose the parameters so 
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Figure 5: Top: Evolution of the length 
of whole bundle including the polar arms 
(the upper curve) and the length of the 
bundle less the polar arms (the lower 
curve). Bottom: The evolution of the 
number of actin filaments sorted out into 
the polar arms. The total number of 
actins, A^, is 1000 in this calculation. 

that Ni ^ Lq holds to assure substantial overlapping of the filaments along the a;-axis. For 
evolution, we define the 'unitary action' by a two-headed myosin: (1) We chose randomly a 
spatial point, say at a; = xm, where the myosin translocates a pair of actins (see (2) below). 
(2) Among all the filaments that extend through the point x = Xm we then choose randomly 
a pair of anti-parallel actin filaments. One of the chosen filament is oriented toward the 
positive X direction (i.e., rightward) and is centered at, say, while the another chosen 
filament is oriented toward the negative x direction (i.e., leftward) and is centered at, say, 
x_. (3) We move these two filaments by the same distance but in the opposite direction, 
according to the scheme described in Fig. ^ As seen from this figure the distance of sliding 
is given by U{x+,x_,xy[) = Min[6, xyi — x+ + £/2, x_ + £/2 — xyi\. (4) As for the rest of 
the filaments in the bundle, those filaments in the region of a; > Xc(x+,a;_) = (x_|_ + xJ)/2 



Figure 4: Snapshots of actin distribution; 
the blue dots indicate the center posi- 
tion of actins oriented rightward, and the 
green ones indicate those oriented left- 
ward. Time proceeds from the bottom to 
top. For the parameter values used, see 
the text. 
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are displaced by +U [—U] if x+ > x_ [x+ < respectively. Also, those filaments in the 
region of x < Xc(x+,x_) are displaced inversely so that U in the last sentence is replaced 
by —U. The evolution of the system is obtained by applying this unitary action from (1) 
to (4) repeatedly. 

The result of the simulation is represented by the distribution of filaments position 
(Figs. ^ and by the density of filaments along x-axis (Fig. |^). The parameters used there 
are = 200, i = 0.25, and Lq = 50. These values are comparable to the experimental 
values (in unit of fim for the lengths) 0. The total time lapse of the evolution is such 
that each filament undergoes, on the average, six times the unitary action of myosin. 
These figures reveal the overall shortening of the assembly of filaments, and also shows the 
clustering of the filaments in rather symmetric fashion with respect to their polarity. In 
Fig. ^ we show another run with extended time lapse. There appear the 'arms' from both 
ends of the bundle, which consist of filaments with unique outward polarities. The length 
of these arms and the number of the filaments in these arms increase in time. Although the 
simulation is restricted in one dimensional space, the characteristic features of the evolution 
thus found are in agreement with experimental observation of (i) shortening of the main 
bundle, (ii) inhomogenization of its thickness and (iii) generation of polar subbundles. 0. 

4 Construction of continuum model 

The above algorithm of simulation can be cast into the evolution equation for the densities 
of actin filaments, p+(x+) and p_(x_), oriented rightward and leftward, respectively. For 
this we needed to assume the smallness of the ratio, b/£ <^ 1, and restrict ourselves to 
the limit of weak spatial heterogeneity, |py <^ £p±. Suppose that a myosin is at x = xm 
and starts to exert the unitary action to a pair of anti-parallel filaments. Let us denote by 
V{x+,X-; xu)dx-^dx- the probability that these two filaments are centered at x = x+ ~ 
x+ + (ix+ (the rightward filament) and at x = x_ ~ x_ + dx- (the leftward filament), 
respectively. In the mean field approximation, this is given as 

^/ p+(x+)p„(x_)e(| - |x+ -XM|)^(f - -xmI) , . 

H[x+,X-,Xm — — n 1— ; — J. -r~, i—r^ TTl — ' y^J 

/|,;_,^l<| dx'+ dx'_p+(xV)p-(x'_) 

where we introduced a step function 6{z) = 1 for 2; > and 6{z) = for z < 0. Us- 
ing this distribution we introduce the weighed average of any function, say C(x+,x_,xm), 
over x_|_ and x_ with xm being fixed, < O >xm5 as < C = J dx+ J (ix_P(x_|_, x_; xm) 
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0{x+,X-,xm)- Then the displacement field u{x;xm) of the actin densities caused by a 
unitary action of myosin at xm is given by u{x,xm) = < U{x+,x^,xm) sgn(x+ — x_) 
sgR{x — Xc{x+,x^))>xf^, where sgn(z) = 29{z) — 1. Using the gradient expansion of the den- 
sities of actin filaments, p+{xm + X^)p_{xm + X_) ~ p+ixM)p^ixM)+ p'+ixM)p~ixM)X^+ 
P^{xm)p'-{xm)X_+ . . . (the prime denotes to take the spatial derivative), the weighed 
average can be evaluated up to the first order approximation to give 



U[X,Xm) 



^ bi f p'4xM, 
£ 6 \p+{xm) P-{xm) J 



sgn(x -Xm)- 



(2) 



Here we have noted that sgn(a; — Xc{x+, can be safely replaced by sgn(a; — xm) in the 
coarse grained description which deals with only the length scales larger than i. The zeroth 
order term —b'^/i in the angular bracket represents the tendency of shortening described 
already, and the second term represents the correction due to the spatial inhomogeneity 
of the filament densities. If, for example, p'^ < and p'_ > hold at x = xu, the latter 
term predicts that the shortening is enhanced compared with the homogeneous case. It 
is understandable since < ^ind p'_ > imply the situations like Fig. |](a) is be more 
likely to be found at xu than those like Fig. |^(b). 

The mean drift velocity of the bundle, v{x), is obtained as the integration of u{x,xm) 
with respect to xu multiplied by the frequency factor k(xm) with which myosins exert 
the unitary actions per unit time and per unit interval along the x— axis. The evolution 
equation for po- (o" = ±) is then ^p^- = — ^ ['^ix)Po-]- Hereafter we the simplest choice 
that the factor k{xm) is an overall constant, say k{xm) = ^o- This case is that one can 
solve the evolution equation analytically and, at the same time, that the essential feature of 
shortening and clustering of the bundle is preserved (see below). Performing the integration 
with respect to xu the evolution equation becomes; 



d_ 

m 



Pa{x,t) 



d_ 

dx 



' ( 2b^ bi^ 
1^0 I — -^x + — log 



P+ix,t) 
p-{x,t) 



Paix,t) 



a 



(3) 



The solution of initial value problem can be given via parametric representation as 
follows: 



x{X,t) =Xe- 



~ H 

66 



1-e- 



log 



p^{x{X,t),t) = p^{X,0) 
Figure |^ shows two examples of the solution of 



dx{X,t) 



dX 



pAx,o) 

P-{X,0) 
a=±. 



(4) 



(5) 



, the one starting from the actin densities 
with in-phase undulation (the left column) and the other one starting with anti-phase 
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undulation (the right column), respectively. If we neglected the logarithmic correction 
term in the solution would simply represent the affine contraction of the bundle, i.e.. 



the correction term acts to promote the clustering of actin filaments of both rightward 
and leftward polarity. Two remarks are in order here: We should note that the simulation 
described in the previous section, as well as the experiments with low myosin concentration, 
would correspond to the slightly different choice of k{xm), that is, k{xm) = i^oLo/L(t), 
where L{t) is the total length of the bundle at time t. This overall factor, L{ty^ would 
change the time scale of the evolution of Fig. |^, but does not change the features of the 
evolution of the density profiles. We would also note that the generation of the arm cannot 
be handled within the present approximation in which the spatial variation of actin densities 
is assumed to be small. 

5 Summary 

We have proposed a simple model for the of contraction of the random uniaxial assembly 
of actin filaments, mediated by the two headed myosin molecules which translocates anti- 
parallel actin filaments. Simulation result have agreed at least in the qualitative level with 
the experimental observation: the shortening of the actin bundle, the clustering of density 
and also the generation of polar arms. The experimental situation studied here may well 
correspond to the stage of evolution where the collective transport by motor proteins had 
first come into existence in the biological world. More generally, it would be interesting 
to study from physicists' viewpoint how a function, in its most primitive form, has been 
first acquired by biological systems at any level of evolutionary history; the problem how 
an allosteric enzyme have acquired the function to translocate the other molecule is a 
challenging problem in this respect. 
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Figure 6: The snapshots of density pro- 
files of the actin filaments oriented right- 
ward (upper curve) and those oriented 
leftward (lower curve), taken from the 
data shown in Fig. |. The time lapse 
is such that the average times of the uni- 
tary actions undergone by each filament 
are, respectively, 0, 1, 2, and 3 from the 
top to the bottom. 




-30 -20 -10 



10 20 30 -30 -20 



10 20 30 



Figure 7: Solutions of the evolution 
equation (^ in the text, from the two ini- 
tial conditions (the top raw). The time 
lapse of the evolved states (the bottom 
raw) are the same for both cases. 
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